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We present the uncertainty relation for the characteristic functions (ChUR) of the quantum me¬ 
chanical position and momentum probability distributions. This inequality is more general than 
the Heisenberg Uncertainty Relation, and is saturated in two extremal cases for wavefunctions de¬ 
scribed by periodic Dirac combs. We further discuss a broad spectrum of applications of the ChUR, 
in particular, we constrain quantum optical measurements involving general detection apertures 
and provide the uncertainty relation that is relevant for Loop Quantum Cosmology. A method to 
measure the characteristic function directly using an auxiliary qubit is also briefly discussed. 
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One might think that everything important has al¬ 
ready been said about the quantum uncertainty of con¬ 
jugate position and momentum variables, discussed for 
the first time almost a century ago m in terms of the 
Heisenberg Uncertainty Relation (HUR). Even though 
the past few years have seen considerable activity de¬ 
voted to describing the uncertainty of non-commuting 
observables in the discrete (mainly in the direction of the 
entropic formulation gHia with emphasis on so-called 
“universal” approach El El HD) or coarse-grained El- 
lie] settings, the most fundamental continuous position- 
momentum scenario appears to be more than well un¬ 
derstood and explored m- For instance, the optimal, 
state-independent entropic counterpart of the Heisenberg 
Uncertainty Relation (HUR) was demonstrated 40 years 
ago |18| , while canonically invariant uncertainty relations 
for higher moments have also been derived m- In atomic 
physics, where the angular momentum of electrons in an 
effective central potential plays a major role, proper mod¬ 
ifications of the uncertainty relation for positions and 
momenta include the relevant eigenvalue of the square 
of the angular momentum operator [20] . and if the elec¬ 
tronic state in question is not the angular momentum 
eigenstate, also the variance of m- In the domain 
of Quantum Electrodynamics the ultimate Heisenberg- 
like uncertainty relations have been obtained for single 
photon states and the coherent states p2l I23| . Even the 
seminal error-disturbance relation by Heisenberg, while 
causing problems in terms of rigorous interpretation, has 
been examined in various different ways |24lE7] . 

Studies devoted to uncertainty relations are often mo¬ 
tivated by a broad network of potential applications. In 
terms of quantum information, for instance, the uncer¬ 
tainty relations have found themselves [55] as important 
ingredients in security proofs of quantum key distribu¬ 
tion |29l ISOj . In experimental studies within the field of 
quantum optics, they have been used in identification of 
quantum correlations such as entanglement |31H33| and 


Einstein-Podolsky-Rosen-steering |34II38| . Beyond quan¬ 
tum information, the uncertainty relations can play an 
important role in various tasks ranging from down-to- 
earth estimation of Hamiltonians |39| to pioneering ex¬ 
periments designed to simultaneously test Quantum Me¬ 
chanics and General Relativity m- 

The present contribution aims to open a new chapter 
in the long history of the uncertainty relations in quan¬ 
tum mechanics. The main subject of our investigation is 
the characteristic function, a notion well known in classi¬ 
cal probability theory. The characteristic function 4> (A) 
related to a probability distribution p (x) is defined as 
the Fourier integral: 

A>{\) = [ dxe^^^p{x). ( 1 ) 

Jr 

The notion of the characteristic function acquires many 
interesting features when considered on the ground of 
quantum mechanics. Let us assume that the probability 
distribution in Q is related to the quantum mechanical 
position space. In this case the characteristic function is 
equal to the average value of the momentum shift 

operator Since this operator is unitary, though non- 
hermitian, it is not an observable. Thus, it might not be a 
natural choice when thinking about uncertainty relations. 
Even when calculating the variance-like quantity 
with A = U — {U'), defined for any unitary operator 
U, one finds that the quadratic term (U"^U) is trivially 
equal to 1. The total uncertainty information is then 
completely contained in the squared modulus |^{7)p. 

Another special feature of $ (A) is related to the mo¬ 
mentum representation. Consider a pure state, so that 
p{x) = \’/{x) / where '/{x) is the position space wave 
function. The momentum wave function ip (p) obtained 
by the Fourier transformation 

iP{p) = ^ fdxe-^P^^^iPix), ( 2 ) 
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naturally provides the probability distribution in momen¬ 
tum space p{p) = |'0 (p) P- Definition ([^ when rewritten 
in the momentum representation gives 

$ (A) = [dptjj* (p) ip{p- hX), (3) 

Jr 

which is the auto-correlation function of the momentum 
wave function. The variable hX, with units of momen¬ 
tum, can be thus interpreted as an auto-correlation pa¬ 
rameter. Obviously, when starting from the momentum 
distribution, an equivalent auto-correlation expression 
can be obtained for the position. To clearly state the dis¬ 
tinction between the two quantum mechanical represen¬ 
tations let us further label by A^, and <i>(Aa;) the argument 
and the function itself, in the case when the characteris¬ 
tic function is calculated for the probability distribution 
p{x). The symbols Xp and $(Ap) shall have the same 
meaning for the momentum density p (p). 

Since the characteristic functions <i>(Aa;) and $(Ap) de¬ 
scribe the auto-correlation of the conjugate variables, and 
considering intuitively that it should not be possible for a 
quantum system to be arbitrarily well localized (or more 
precisely, well auto-correlated) in both position and mo¬ 
mentum variables, one would expect that there must ex¬ 
ist some constraint on these functions, in the form of a 
quantum mechanical uncertainty relation. 

By construction, the modulus of the characteristic 
function is trivially upper-bounded by 1. In general, 
however, it is hard to provide a non-trivial, A-dependent 
upper bound. In the mathematical literature one can 
find results relying on the assumptions that the prob¬ 
ability distribution has finite support, a given variance 
or is bounded im na. More elaborate studies take 
into account higher moments or even the entropy ga¬ 
lls] (other relevant mathematical references can be found 
in gSlIlTI). A single restriction does not however lead to 
an upper bound that is sharper than 1 , so one needs to 
assume more gj. Unfortunately, both the finite support 
and the upper bounded value of p{x) do not work well 
even for the most basic quantum-mechanical wave pack¬ 
ets, such as the family of Gaussians. For example, the 
variance-dependent upper bounds (with unbounded sup¬ 
port of the distribution) given in Ref. gT] are of the form 
|$(Aa;)| < with g{ax) = CX^/ |Aa;| -I-tt), 

where (Ta, is the standard deviation, is the maximum 
of the distribution and C is a numerical constant. In 
principle, this upper-bound might be useful in the cur¬ 
rent context, provided that one is able to apply also the 
HUR and bound ax > h/{2ap). However, it is clear that 
the HUR cannot be utilized since the upper bound is a 
monotonically increasing function of ax , and does not im¬ 
ply |$(Aa;)| < In addition, the optimized 

(with respect to ax and Up) sum of and its 

momentum counterpart gives only a trivial bound. 

An easier scenario occurs when dealing with lower 
bounds for the characteristic function, due to a well 





FIG. 1: (Color online). The upper bound I?(?iAa;Ap) (solid, red 
line) with its upper (dotted, equal to 2) and lower (dashed, 
equal to 1) values emphasized. The blue dashed-dotted line 
represents the left hand side of ([^ calculated for the Gaussian 
state discussed after Eq. Q. 


known fact (cf. for example Lemma 1 of |41|1: 

|$(A)| >Re$(A) > 1-^AV^ (4) 

The above inequality holds for any A, provided that the 
variance a^ is finite. However also in this case, the HUR 
cannot be directly applied as it would increase the lower 
bound in Q . Let us mention that the characteristic func¬ 
tion has been considered in gS] , mainly in the context of 
time-energy uncertainty relations along the lines of the 
relation 0 . 

In the current discussion we aim to go beyond moment 
approximations (and related bounds) and take into ac¬ 
count the complete information content stored in 
We thus present the uncertainty relation for the position 
and momentum characteristic functions (ChUR): 

Theorem 1. The sum of squared moduli of the po¬ 
sition and momentum characteristic functions is upper 
bounded: 


mxx)\^ + mxp)\^ <B{hXxXp), ( 5 ) 


by 


B{-i) = 2V2 


y/2- cos ( 7 ) 
1 -I- cos ( 7 ) 


( 6 ) 


The characteristic function is a linear functional of the 
probability distribution, so its modulus squared is con¬ 
vex. Since any quantum state can be written as a convex 
sum of pure states, the ChUR (|^ is satisfied by all mixed 
states. Moreover, the above result can be immediately 
generalized to the multidimensional scenario in which Q 
is defined as a d^x integral with e^^'^. In that case, Eq. 
[^remains valid with A^, (Ap) relpaced by (Ap), and 
consequently Aa;Ap replaced by the scalar product Xx- Xp. 

We will provide an outline of the proof below. First, 
let us discuss the function B {■) in inequality ([^, which 
depends on the dimensionless parameter hXxXp, and is 
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plotted in Fig. The upper bound is clearly peri¬ 
odic and varies between 1 and 2 , the latter being the 
trivial value. When hX^Xp = 2kTT for any integer k 
we have B{2k-K) = 2, and Eq. (§ gives no restriction 
on the sum of characteristic functions. This fact is a 
non-trivial emanation |49| of the full commutativity of 
spectral position and momentum projections |5()| . In 
other words, the trivial bound equal to 2 is saturated 
by a distribution p [x] that is a (normalized) version of 
the Dirac comb (wave function being the proper limit of 
the sum of shifted Gaussians) with period 271 jXx- The 
characteristic function $(Ax) is equal to 1 in this case. 
At the same time, $(Ap) is the auto-correlation func¬ 
tion of the comb with the related correlation parameter 
hXp. In this case hXp = 2 A: 7 r/Aa;, so this parameter fits 
the comb period, giving $(Ap) = 1. In the “opposite” 
situation, when hX^Xp = {2k -I- 1 ) tt the upper bound 
reaches its minimal value 1. This case captures the maxi¬ 
mal non-compatibility of the position/momentum couple 
since both wave functions i/' (x) and ijj {p) cannot be si¬ 
multaneously well auto-correlated. 

Before we prove Theorem 1 we would also like to show 
that the UR © is stronger than the HUR. To this end 
we parametrize A^, = ^fajh and Xp = y/ab/h, where a is 
a non-negative dimensionless parameter and b denotes a 
parameter with unit of position. By taking the square of 
inequality (|^ and neglecting the positive contributions of 
order cr^, one immediately arrives at |<i)(Aa;)p > 1 —A^cr^, 
and the similarly for |$(Ap)p. This lower bound together 
with the above parametrization weakens inequality (|^ to 
the form 

2-a (b-^al + {b/hf cr| j < S(a). (7) 

Since for a > 0, B{a) = 2 — a + 0{a?), the above relation 
divided by a implies in the limit a —>■ 0 that 

\<b-^al + {b/hfal. ( 8 ) 

The minimum of the right hand side occurs for b = 
y/haxIcTp and gives the HUR. Thus, the ChUR (|^ is 
strictly stronger than the Heisenberg Uncertainty Rela¬ 
tion. By a straightforward calculation one can check that 
the left hand side of in the case of ^^{x) being a Gaus¬ 
sian state and with all the above assignments (for Xx, Xp 
and b) is equal to 2e““/^. In the limit a —)■ 0, the Gaus¬ 
sians saturate the bound as can be seen in Fig. [l] 

Proof. Let us now outline the proof of Theorem 1. We 
start by taking three vectors |^i) = |t|/), [^ 2 ) = |4>) 

and 1 ^ 3 ) = |'I') with |'I') being normalized. The 

positive semi-definite, hermitian Gram matrix of this set 
of vectors is equal to 


/ 1 <i>(A,) ^Xp)\ 

G= $*(A,) 1 U , 

V l>*(Ap) n* 1 J 


where U = (v]>| |4>). The condition of positive- 

semi-definiteness of G leads to a single nontrivial inequal¬ 
ity det G > 0 which explicitly reads 

l-A-|U|V0-f0* >0, (10) 

where (we omit here the A arguments) A = -|- 

and 0 = Gonsider now the parity transformation 

Xx >—t —Xx and Xp 1 -^ —Xp. A basic property of the char¬ 
acteristic function is that $(—A) = $*(A). Moreover, the 
well known Baker-Gampbell-Hausdorff formula (equiva¬ 
lent to the Weyl commutation relations) 




( 11 ) 


provides the transformation rule U 1 —>■ e xhe 

terms A and |U|^ present in (10 1 are thus invariant 


with respect to the above t rans formation, while 0 1 —>■ 
g-thx^XpQ* ^ Gbviously, Eq. (10 1 must also hold for the 
transformed quantities. 


If we now take the arithmetic mean of (10 1 together 


with its transformed counterpart, we will obtain an in¬ 


equality of almost exactly the same form as ( 10 ), with 


the only difference that 0 is now multiplied by a com¬ 
plex constant Z (i.e. 0 1 —>■ Z0 ) of the form Z = 

|(1 -I- In the final steps we resort to the fact 

that Z 0 < |Z| | 0 | and the arithmetic-geometric mean 
inequality |0| < 1 |U| A. The same procedure is applied 
to the second, conjugated term Z*0*. The above deriva¬ 
tion leads to an inequality 1 — A — |U|^ -|- \Z\ |U| A > 0 
which after a single rearrangement is brought to the form 


A < 


1-I»l^ 

1-|Z||U|- 


( 12 ) 


Since the parameter 0 < |r2| < 1 can in principle assume 
any value, we maximize the right hand side of the above 
inequality with respect to it. The global maximum is 

found at |U| = (1 — 1 — \Z\^)/ jZj, and leads to the 

final result presented in Theorem 1, where the identity 
\Z^ = i[l -|- cos{hXxXp)] has been utilized. □ 


Let us mention that if one starts the proof with an al¬ 
ternative choice 1 ^ 2 ) = ^ I'l') 1 ^ 3 ) = p I'k), the coun¬ 

terpart of inequality ( [To| is equivalent to the Robertson- 
Schrodinger uncertainty relation. Since in the current 
case we deal with non-hermitian displacement operators, 
such trivial correspondence does not occur. 

ChUR and Quantum Optics — Though the displace¬ 
ment operator does not correspond to an observable, its 
mean value can be measured directly in a quantum op¬ 
tical experiment if an ancillary qubit is used. Gonsider 
the quantum circuit shown in Fig. The qubit is ini¬ 
tialized in the |-|-) state, where |±) = (| 0 ) ± \l))/x/2 are 
the eigenstates of the Pauli operator ax. The quantum 
system of interest is in the arbitrary state p. The logic 
gate is a controlled displacement operator, defined by 

|0)(0|i + |l)(l|e*^-’^ (13) 


( 9 ) 
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FIG. 2: Quantum circuit to measure the characteristic func¬ 
tion of p{x) directly. 


where the displacement is along the x direction. Detect¬ 
ing the qubit in the ± basis gives the probabilities 

P± = ^(1 ± (cos(App))p). (14) 

Detecting the qubit in the dj, basis, defined by eigenstates 
|±z) = (|0) ± i\l))/y/2, gives 


By calculating the ordinary Fourier transform [like in 
Eq. ([^, but with h = 1 and the conjugate param¬ 
eter A in units of inverted position] of both detection 
probability functions we obtain Q(A) = M(A)$(A) and 
'P(A) = M(A)$(kA). We are thus in position to propose 
a general uncertainty relation for detection masks: 

[dy (|Q(y)P + \V{y)\^) < [dX\M{X)\^B (hnX^) . 

Jm. Jm. 

(18) 

The above result should prove useful in studies devoted 
to quantum aspects of EPR-based ghost imaging m and 
security protocols for compressive quantum imaging |55l 
l56] . A "periodic" variant of Eq. ( [I^ (when M{x) is 
a periodic function) is already applied in experimental 
entanglement detection with periodic amplitude masks 


P±i = T (sin(App))p). (15) 

Then the characteristic function can be obtained by com¬ 
bining the measurement results, since (e®'''p^)p = P+ — 
P_ — iP-i + iP+i- This type of qubit-assisted measure¬ 
ment scheme can be realized using existing technologies 
in several different platforms Eni^. 

As an application of Theorem 1 in Quantum Optics 
we describe the mutual incompatibility of measurements 
made in position and momentum space using detectors 
with arbitrary apertures. In our analysis we discuss a 
very general model, in which the detection aperture is 
described by an arbitrary transmittance function M (x). 
For single photons, for example, this detection scheme is 
implemented by the propagation through a general am¬ 
plitude (and phase) spatial mask modelled by A{x) = 
A{x)e‘^^^\ followed by its subsequent measurement with 
a full multi-mode detector |53| . The aperture function 
0 ^ A{x) ^ 1 provides M{x) = \A{x)\‘^ = A‘^{x) while 
the mask phase profile (j>{x) does not affect the transmit¬ 
tance. The probability Q (y) that the quantum particle 
is detected with this mask function is then given by 

Q{y) = [ dxM{x + y)p{x), (16) 

Jm 


where y can be thought of as the location parameter that 
defines the mask. For example, if the mask is an aperture 
of size d (in some experimentally relevant units), then 
M (x) might be chosen to be equal 1 for 0 < x < S and 0 
elsewhere. In that case (161 is simply the probability of 
finding a quantum particle on the interval [—y, 6 ~ y]- 
A similar construction can be done in the momentum 
picture. One only needs a parameter k mapping the mo¬ 
mentum variable into the position space, so that Kp is a 
position-like variable. To have the readout function de¬ 
pending on the position-like variable we define the coun¬ 
terpart of (161 as follows 


The derivation of (181 is very simple. Due to Parse- 
val’s theorem from Fourier analysis, the left hand side 
translates directly to the A-domain. To obtain the right 
hand side we apply the UR for characteristic functions 
from Theorem 1, with X^ = X and Ap = kX. It is worth 
mentioning that Eq. remains valid for any complex¬ 
valued function M{x). 

ChUR and various theories in physics — Besides its 
fundamental interest, the ChUR derived in this paper is 
related to several issues across various fields of physics. 
First of all, our approach remains valid with $(Aa;) and 
$(Ap) substituted by ('I'| U I'k) and (ikl W |'k), whenever 
the unitary matrices U and W satisfy the Weyl-type com¬ 
mutation relations UW = e^^WU. A prominent example 
provided by Schwinger m and given hy (j) = 2 tt jd, where 
d is the dimension of the Hilbert space, is a sort of prereq¬ 
uisite for the fruitful theory of Mutually Unbiased Bases. 

Moreover, since Theorem 1 involves operators of the 
form it becomes valuable when the operator O does 
not exist itself (consequences of so called Stone-von Neu¬ 
mann theorem). A particularly interesting example of 
the number-phase uncertainty 1581150] (phase operators 
are not well defined) has just been described [OT] along 
the lines of Theorem 1. Here we would like to briefly 
touch upon the broad theory of Loop Quantum Grav¬ 
ity [T^, in which the so called Ashtekar connection jOO] 
m plays the role of the canonical "position" vari¬ 
able in a field-theoretical sense m- While moving to a 
quantum description, the problem appears as there is no 
local operator A\^{x), and one needs to resort to unitary 
holonomies. The standard approach to quantum uncer¬ 
tainty relations cannot thus be directly applied, however 
one can involve the Weyl algebra [M] and use Theorem 1. 

To explain better the idea behind the above prescrip¬ 
tion, we would like to discuss a very simple case from the 
field of Loop Quantum Cosmology. To this end we start 
with the well known FLRW metric 


Piv) 


dpM [up + y) p{p). 


(17) 


= —(Adt^ + a^(t)dS^, (19) 
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where a{t) denotes the dimensionless scale factor and S 
refers to the 3-dimensional space. We further recall two 
time-dependent variables |65]: b = dja and V = Vga^, 
denoting the Hubble parameter and the physical volume 
of the expanding Universe respectively (Vb is the coordi¬ 
nate volume). These variables satisfy the following Pois¬ 
son bracket relation m 

{b, V} = ±4 ttG/c'^ = Q, (20) 

where the ± sign depends on the orientation and is ir¬ 
relevant in our considerations. If the operator b existed, 
then ( [20| would lead us to the UR: abcry > iirfiG/c^ as 
stated~j^ in Eq. 11.21 of |^. Since the Ashtekar con¬ 
nection of this well-studied model is given by A\ = /3a(5* 
177], it becomes obvious that b cannot be promoted to a 
quantum mechanical operator. As this limitation is not 
shared by the holonomy Ub{Xb) = one can use Ub 

together with and apply Theorem 1. In particular, 

if we set Xy = "^/{fiQXb), so that the bound B is equal 
to 1, and use Q to extract the variance Uy, the ChUR 
provides the uncertainty relation of the form: 

^Xb\{Ub{Xb))\<av. ( 21 ) 

The above UR is a formally right way of bounding the 
fluctuations of the volume of the Universe in terms of the 
volume shift operator Ub, relevant for understanding of 
the big-bang singularity j^. Note also that this exam¬ 
ple actually represents Quantum Mechanics subject to 
the Bohr compactiflcation. In other words. Theorem 1 
is the only path towards URs in theories (such as Loop 
Quantum Cosmology |67j) with the Bohr compactiflca¬ 
tion involved. 

Looking into future, a further development of the dis¬ 
crete counterpart of the presented theory might bring 
useful results, for instance, in compressed sensing, as the 
Dirac comb state has no counterpart in various discrete 
systems. Generalizations of URs for the electromagnetic 
held (as discussed in |68|) might bring a better physical 
insight into the role played by Gauss linking numbers, 
or contribute to a better understanding of quantum ef¬ 
fects for the gravitational held in a hot universe |69| . 
We also believe that our approach will be influential to 
the theory of quantum optical characteristic functions. 
Questions about non-classicality of light are being asked 
and studied m in terms of the characteristic P-function 
reconstructed from the data accessible in experiments 
ini, with a relevant Altering procedure based on auto¬ 
correlations EH- 
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